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We derive analytic solutions of a chameleon scalar field <j> that couples to a non-relativistic matter 
in the weak gravitational background of a spherically symmetric body, paying particular attention 
to a field mass vtla inside of the body. The standard thin-shell field profile is recovered by taking 
the limit mAr c — * oo, where r c is a radius of the body. We show the existence of "no-shell" 
' solutions where the field is nearly frozen in the whole interior of the body, which does not necessarily 

, correspond to the "zero-shell" limit of thin-shell solutions. In the no-shell case, under the condition 

■ mAr c 3> 1, the effective coupling of cj> with matter takes the same asymptotic form as that in the 

thin-shell case. We study experimental bounds coming from the violation of equivalence principle 
as well as solar-system tests for a number of models including f(R) gravity and find that the field 
is in either the thin-shell or the no-shell regime under such constraints, depending on the shape of 
' scalar-field potentials. We also show that, for the consistency with local gravity constraints, the 

field at the center of the body needs to be extremely close to the value 4>a at the extremum of an 
effective potential induced by the matter coupling. 
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I. INTRODUCTION 



Recently there have been a lot of efforts to understand the origin of dark energy under the framework of modified 
gravity theories (see Refs. [l| for reviews). Presumably the simplest modified gravity model of dark energy is so-called 
^ ' f{R) gravity in which / is a function in terms of a Ricci scalar R (2j. If we transform to the Einstein frame by a 
conformal transformation, f(R) gravity is described by a scalar-field action with a potential that has gravitational 
origin Q. In the Einstein frame, the scalar field couples to a non-relativistic matter with a coupling of the order 
of unity (Q = — 1/V&) [4]. For such large coupling models the potential generally needs to be carefully designed 
to satisfy both cosmological and local gravity constraints. See Refs. [!, 0, 0, § for the construction of viable f(R) 
OP | models. 

Scalar-tensor models of dark energy [T3] also give rise to a similar coupling Q in the Einstein frame, where the 
strength of Q depends on the gravitational coupling with the scalar field in the Jordan frame [lp. Thus, in modified 
gravity models of dark energy, it is crucially important to appropriately study the compatibility of such couplings 
with local gravity experiments as well as with a late-time acceleration of the Universe preceded by a standard matter 
era. 

It is also known that there are many massless scalar fields such as a dilaton in superstring and supergravity theories 
[IH ■ In general, these fields couple to matter as well as gravity with strengths of the order of unity. The existence of 
such scalar fields gives rise to a strong violation of equivalence principle, thus posing a difficulty for the compatibility 
with local gravity experiments in the absence of field potentials. One of the attempts to tackle this problem is the 
so-called run-away dilaton scenario [l3| in which the field gradually decouples from matter in the large-field regime 
under an appropriate choice of the dilaton coupling. It is then possible to apply the run-away dilaton scenario to dark 
energy provided that the dilaton has a monotonically decreasing potential toward larger <fi [141 ]. 

Another very interesting attempt to reconcile large coupling models with local grav ity constraints is to use 
"chameleon" scalar fields whose masses depend on the environment they are in (Til. Il6j. The chameleon mecha- 
nism naturally works for monotonically decreasing potentials that appear in the context of dark energy without 
adjusting the strength of the matter coupling Q [171 ] . The chameleon mechanism has been applied to many models- 
such as self- interacting scalar field models [Til], f(R) gravity frj| . JT9I [2fj| . scalar-tensor theories 0, HT|, [13] , curvature 
couplings to a scalar field (23| , radions in braneworld models |24|. and theories whose matter coupling is much larger 
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than the gravitational coupling (2f|. The presence of the matter coupling induces an extremum of the field potential 
around which the field can stay. If the density of the matter is sufficiently high as the interior of a compact object, the 
field acquires a heavy mass about the potential minimum so that the violation of equivalence principle is suppressed 
even if Q is of the order of unity. Meanwhile the field has a lighter mass in a low-density cosmological environment 
relevant to dark energy so that it can propagate freely. 
In Refs. [H [3 

it was shown that local gravity constraints can be satisfied when a spherically symmetric body has 
a thin-shell. When the thin-shell is formed the chameleon field stays around the potential extremum at 4> — 4>a inside 
of the body in the region < r < 7*1 , where n is close to a radius r c of the body. The field evolves in the thin-shell 
region characterized by r\ < r < r c due to the dominance of the matter coupling term Qpa, where pa is a mean 
density of the body. As long as the condition, (r c — r\)jr c -C 1, is satisfied, the effective coupling Q c g outside of the 
body becomes much smaller than the order of Q so that the models are consistent with local gravity experiments. 
See Refs. [2r| for recent experimental constraints on model parameters under the chameleon mechanism. 

The analysis in Refs. [HlQil assumes that the field is frozen (0 = 4>a) in the regime < r < r\ without explicitly 
taking into account the field mass tua- In this paper we derive field solutions in three regimes (a) < r < n, (b) 
n < r < r c , (c) r > r c , and determine integration constants analytically by imposing appropriate boundary conditions 
at r = 0, 7*1, r c , 00, respectively. Under the condition that the field mass tub about the potential extremum far outside 
of the body is nearly massless, the field profile takes fairly simple forms as given in Eqs. (|35|) - (|37|) . 

We will show that standard thin-shell solutions found in Refs. [HI, [f6| correspond to the limit mAr c — > 00 with 
(r c — ri)/r c <C I. There exist "no-shell" solutions where the effect of the Qpa term does not become important in 
the whole interior of the body. These are wider class of solutions than the "zero-shell" limit of thin-shell solutions as 
we show later. The effective coupling Q cff outside of the body in the no-shell case takes the same asymptotic form 
as that in the thin-shell case in the limit rriAT c — > 00. Hence it is possible to satisfy local gravity constraints even if 
the body does not have a thin-shell. We also study experimental bounds for concrete models and show that the field 
profile consistent with these bounds corresponds to either thin-shell or no-shell solutions with uiat c ^> 1 , depending 
upon the form of field potentials. 

This paper is organized as follows. In Sec. |TT] we provide a number of concrete models which we will use in a 
later section. These are basically motivated by the potentials that appear in the context of dark energy. In Sec. IIIII 
we review the chameleon mechanism and derive new solutions by properly considering the field profile in the region 
< r < 7*1. In Sec. HVI we obtain three different solutions: (i) thick-shells, (ii) thin-shells and (iii) no-shells, by taking 
some limit for the solutions derived in Sec. IIIII In Sec. [V]we apply our formula to a number of field potentials and 
clarify the regime of the field profile when local gravity constraints are satisfied. Sec. I VII is devoted to conclusions. 



II. MODELS 



We consider a scenario in which a scalar field <f> with potential V((j>) couples to a non-relativistic matter with a 
Lagrangian density C m . The action we study is given by 



S = / d A x*J=g~ 



Ml 



I 



-^R--(y4>Y~V{4>) 



(1) 



where g is a determinant of the metric g^ u , M p \ = I/VsttG is a reduced Planck mass (G is gravitational constant), R 
is a Ricci scalar, and are matter fields that couple to a metric g)fl related with the Einstein frame metric g^ u via 



7 (0 - p2Q, 



9tiv ■ 



(2) 



Here Qi are the strength of couplings for each matter field. In the following we shall use the unit M p \ = 1, but we 
restore M p \ or G when it is needed. 

There are theories that give rise to a constant coupling Q for each matter field. Let us consider the following 
scalar-tensor action [T(| 



S 



d x^J-g 



\f{$)R - \{\ - &Q 2 )F (0)(V0) 2 - U(4>) 



d 4 x£ m (^ m ,g^) , 



(3) 



where a tilde represents quantities in the Jordan frame and 
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Under a conformal transformation. 



9^ = F(<l>)g»v , (5) 
we get the Einstein frame action ^ with the field potential 

1/(0) = U^/F^f . (6) 
It is worth mentioning that the action ([3]) is equivalent to the following Brans-Dicke theory (2?| with a potential U((f>): 



S 



(Vx) 2 -PWx)) 



(7) 



where x = F{4>) = e 2< ^ an d wbd is a Brans-Dicke parameter that is related with the coupling Q via the relation 
3 + 2cjbd = 1/(2Q 2 ) The /(i?) gravity corresponds to the coupling Q = -1/V&, i.e., wbd = [H. 

If the field potential is absent, solar-system experiments constrain the Brans-Dicke parameter to be wbd > 4.0 x 10 4 
[29l ]. This translates into the bound: \Q\ < 2.5 x 10~ 3 . When \Q\ is of the order of unity, as in the case of f(R) 
gravity, the presence of the potential is essentially important for the consistency with local gravity experiments. If 
the scalar field is responsible for dark energy, it is natural to design the potential V((f>) so that the field runs away 
toward larger \(f>\ with a tracking behavior at late times. An example of this type is the inverse power-law potential 

M 



V(4>) = M 4+n (/)- 



(8) 



where M has a unit of mass and n is a constant. The scalar field is almost trapped at the extremum of an effective 
potential induced by the presence of the matter coupling Q with \<f>\ much smaller than M p \ (as we will see later). As 
long as the field stays at the extremum in the region of high density, it is possible to satisfy local gravity constraints 
through the chameleon mechanism. 

In the context of f(R) gravity, Hu and Sawicki Q and Starobinsky 0] proposed models that can be consistent with 
cosmological and local gravity constraints. In the region where R is much larger than a critical value R c (R c is of the 
order of the present cosmological Ricci scalar), these models correspond to the Jordan frame action 



with 



S = J d 4 aVH^/(£) - / d 4 x£ m (* rn ,^), 



f(R) = R-fiR c [l-(R/R c )- 2n ] 



(9) 



(10) 



where \x and n are positive constants @. Comparing the action @ with ^ we find that the potential in f(R) gravity 
(Q = — l/v6) corresponds to U = (FR — f )/2 with a dynamical scalar field <f> = (\/6/2) InF, where F is related with 
/ via the relation F = df jdR [TjJ. The field potential for the model (|10[) is given by 



U(<f>) 



yR c 
2 



2n 



(11) 



For the action ([3]) with arbitrary couplings Q, one can also construct viable models by generalizing the potential 
(fTTj) in f(R) gravity. An explicit example of this type of potential is given by [TJj] 



U(cf>) = V [1 - C(l - e- 2 ^f] , (12) 

where Vq > 0, C > 0, < p < 1. The model (fTO)) in f(R) gravity corresponds to p = 2n/(2n + 1). The scalar field 
mass gets larger as p becomes closer to 1. Since \(f>\ is much smaller than 1 in the region of high density, the potential 
V(4>) = e 4< ^^U(<p) in the Einstein frame is almost identical to the potential U(4>) in the Jordan frame for \Q\ < 1. 

If the action JT]) does not originate from the scalar-tensor action, it can happen that the strength of the coupling 
Q is not the same between different matter species. For example, in coupled dark energy scenario in Ref. [ll|, it is 
assumed that a quintessence field couples to dark matter but not to baryons. 
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III. CHAMELEON MECHANISM 



In this section we revisit the chameleon mechanism paying particular attention to the scalar-field mass inside of a 
spherically symmetric body. The contribution of metric perturbations to the scalar-field equation is neglected as we 
only consider a weak gravitational background. As a matter source we take into account a non-relativistic fluid whose 
pressure is negligible relative to its energy density. 

Varying the Einstein frame action ([T]) with respect to the field (j>, we get 

i 

where Tf$ = (2/ \/—g^)5C m /8g^ u is the energy momentum tensor of the i-th matter. The trace of the i-th matter 
is given by TM = g^T$ = —pi for a non-relativistic fluid, where pi is an energy density. 

From the comparison of Eq. j2]) with Eq. ([5]), it is clear that /?, has a meaning of the energy density in the Jordan 
frame if the action ([T]) originates from ^ under the conformal transformation. The energy density in the Einstein 

frame corresponds to p\ E ^ = pte 4 ®^, but this does not satisfy the usual continuity equation. It is more convenient 
to introduce the quantity pi = pi e 3 *^, which is conserved in the Einstein frame 1 . Then Eq. Q13p reduces to 

i 

In a spherically symmetric background Eq. yields 

fl + 2d^^dV^_ (15) 
dr 2 r dr d(j> ' 

where r is a distance from the center of symmetry and the effective potential V c r is defined by 

v oS (<p) = v^) + Y,p^ ■ ( 16 ) 

i 

In the following, unless otherwise stated, we shall consider the case in which couplings Qi are the same for each matter 
component, i.e., Qi = Q and pi = p. 

When Q > the potential has a minimum for the models with < 0. For example, the potential ([5]) gives rise to 
a minimum at 



4*M 



n f M ^ 4+ " 



P \M pl 



i/(,i+i) 

Afpj , (17) 



where we recovered the reduced Planck mass M p \. The condition, <f>nj -C M p \, needs to be satisfied for the consistency 
with local gravity constraints [3, [l6[ ■ 

When Q < 0, as in f(R) gravity, the potential has a minimum for the models with V.^, > 0. In fact the potential 
given in Eq. Q2J), i.e., V(<j>) = Vbe 4C ^ [l - (7(1 - e~ 2 Q<t>)P], satisfies the condition vj > for Q < 0. Then the 
effective potential V c ff has a minimum at 



1 (2pCV^ 1,(1 ~ p) 



which exists in the region <pM < 0. Note that the order of Vq is not much different from the present cosmological 
density po ~ 10 _29 g/cm 3 if this potential is responsible for the accelerated expansion today. Hence in the region 
p 3> Vq the condition \4>m \ «C M v \ is well satisfied [Io| . 

In the following we assume that the spherically symmetric body has a homogeneous density p = pa and that the 
density is homogeneous with a value p = ps outside of the body. The mass of this body is given by M c = (47t/3)/?a?~c; 



1 In the Friedmann-Robertson- Walker cosmological background this means that pi satisfies the equation pi + 3Hpi = (H is a Hubble 
parameter), while the equation for p' E ' is p;' E ' + 3Hp^ = Q;0p' B ' f3ltl . 
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where r c is a radius of the body. The effective potential V e $ has minima at field values <pA and 4>b characterized by 
the conditions 

V40 A ) + Qp A e Q * A = 0, (19) 
V<p(^ B )+QpBe Q * B =0. (20) 

The former corresponds to the region with a high density (interior of the body) that gives rise to a heavy mass squared 
m A = d d 02 ff {4>a)i whereas the latter to the lower density region (exterior of the body) with a lighter mass squared 

We impose the following boundary conditions: 

^(r = 0) = 0, 4>(r -> oo) = 4>b ■ (21) 
dr 

Equation (| 1 5j) shows that we need to consider the potential (— V e s) in order to find the "dynamics" of <j> with respect 
to r. This means that the effective potential (— V c f-r) has a maximum at <p = 4>a- The field <j) is at rest at r — and 
begins to roll down the potential when the matter-coupling term QpAe®^ becomes important at a radius r\. 

If the field value at r = is close to <\>a, the field stays around <j)A in the region < r < r\. The body has a 
thin-shell if r\ is close to the radius r c of the body. Whether >'i can be dose to r c or not also depends on the negative 
mass squared (—m A ) at the potential maximum. If tua is much larger than l/r c , it is likely that we need to determine 
the boundary condition of 4> at r = very close to 4>a in order to obtain thin-shell solutions. 

In the region < r < ri, the r.h.s. of Eq. (|15p can be approximated as dV e f{/d(j> ~ ra 2 A {4> — 4>a) around 4> = 4>a- 
Then the solution of Eq. (|15p is given by <j>(r) = <Pa + Ae~ mAT /r + Be mAr /r, where A and B are integration constants. 
To avoid the divergence of <p at r = 0, we require that B = —A. Then the solution is 

A(e~" lAT — e" 1Ar ) 

0(r) = <f>A + — (0<r<n). (22) 

r 

This automatically satisfies the boundary condition: 3^(r = 0) = 0. 

In the region r\ < r < r c the field \4>(r)\ evolves toward larger values with the increase of r. Since \V^\ -C \Qpa&^^\ 
in this regime one has dV c f[/d<p ~ Qpa in Eq. (|15p . where we used the condition -C 1. Hence we obtain the 
following solution 

<t>(r) = \Qpat 2 --+D ( ri <r<r c ), (23) 
6 r 

where C and D are constants. 

In the region outside of the body the field \<f>\ climbs up the potential hill toward larger values. The kinetic energy 
of the field dominates over the potential energy, which means that the r.h.s. of Eq. (TT5)) can be neglected relative 
to each component of the l.h.s. of it. Taking into account the mass term dV e ^/d(j} — rn1(<p — 4>b) on the r.h.s. of 
Eq. (|15[) . we obtain the solution <f>(r) = 4>b + Ee~ mB ( r ~ rc ' > jr + F e m B(r-r c ) j r integration constants E and F. 
Demanding the boundary condition, cj>(r — » oo) = <fig, one has F = and hence the solution is given by 

»-ms(l — r a ) 

4>{r) = 4> B + E (r>r c ). (24) 

r 

We match three solutions (|22p . ([2^)1 and by imposing continuous conditions for (j> and d(f>/dr at r = r\ and 
r = r c . Then four coefficients A, C, D and E are determined accordingly: 

c = gi£2[(0g - <M + (rf - rg)Q PA /6] + [agrfCg^^ - e^) - g ir c 2 ]Qp A /3 

m j4 (e- mAri + e m - 4ri ) s 2 - m B si ' ^ ' 

^ = --(C + Qp A rl/3), (26) 
si 

£ = - — (C + Q PA rl/X), (27) 

£> = 4>B-Qp A r 2 J&+-{C + E), (28) 
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where 



s a = m Ari (e- mAri +e mAri ) + e 
s 2 = 1 + m B r c . 



(29) 
(30) 



Especially when the conditions, msT c <C 1 and m A 3> mj, are satisfied so that the contribution of the ma- 
dependent terms are negligible, we have 



A = — 



1 



m A (e- mAri + e mATl ) 



1 



?b - <I>a + ^Qpa{t\ - r 2 c " 



C = n 



-m A r 1 _ p m A r 1 



1 



m A r\{e- mAri + e m ^) 



1 



D = <Pb - -jQpAr, 



2 

c ' 



E = -n( 



'B - <Pa 



) - -Qp A ii ( 2 



ri 



-m A ri _ .mjri 



TOA(e 



</>B - <f>A + -^QPA{r\ - f^) 



(31) 

(32) 

(33) 
.(34) 



This leads to the following solution 



1 



m A (e- mAri + e m ^''i) 



5B - 4>A + ^QPAirx - r c 



-m A r _ „m A r 



(0 < r < n), (35) 



(r) = 0b + ^PA(r 2 - 3r c 2 ) + 9£gL 



-m A rx _ „m A r\ 



1 



b B -cl> A + -Qp A (rl-rl) 



r i i \ 
— (n <r < r c ), 

r 



(36) 



<f>(r) = <p B - 



6 



ri{cj)B - <i> A ) + -^QpArl [ 2 + ^ ) I 1 - - 



ri 



!>b - 0a + -^Qpa{tI - r 2 c ] 



-m. B (r-r a ) 



(r > r c ) . (37) 



While a similar matching procedure given above was considered in Ref. [31(, the analytic expression (|3"5|) - (f3"T)) has not 
been derived so far. 

The radius r\ is determined by the following condition 



m A W r l) - <Pa] = QPA, 



which translates into 



9B - <\>A + T,QPA{r\ - r\) = 



6Q$ C m A r 1 {e mAri + e~ mAri ) 
{m A r c ) 2 



a rn A ri p — m A ri 



(38) 



(39) 



where $ c = AI c /(8irr c ) = p A r^/6 is a gravitational potential at the surface of the body. Under this relation the field 
profile (f3"T)) outside of the body can be written as 



(r) = 4> B - 2Q$ c r c 



rl „n 1 ( mjinie" 1 ^ 1 + e~ mAri ) 



1-4 + 3 



r'c r c (m A r c ) Q 



= -rn B (r-r e ) 



(r > rc) ■ (40) 



IV. THICK-SHELL, THIN-SHELL, AND NO-SHELL SOLUTIONS 



In this section we derive three different types of solutions: (i) thick-shells, (ii) thin-shells, and (iii) no-shells by 
taking appropriate limits for the results derived in the previous section. 



A. Thick-shell solutions (n — > 0) 



If the field value at r = is away from 4>a, it happens that the field rolls down the potential for r > 0. This is the 
thick-shell regime in which the field inside of the body (0 < r < r c ) changes significantly. 
More precisely the following relation holds in the thick-shell case: 

m\ [0(0) - <f> A ] > Qpa ■ (41) 

The field value at r = is given by 0(0) = 0b — (l/2)QpArl by taking the limit r\ — > in Eq. (|35j) . Then the 
condition (|4Tj) can be expressed as 

£th> ^ + Kkr (42) 

where eth is the so-called thin-shell parameter [l~5|, [l6| defined by 

_ <t>B ~ <j>A 

eth -^g$T- (43) 

Equation shows that, in the thick-shell case, eth is larger than the order of unity. 
Taking the limit n -> in Eqs. and l[3"7|). we get 

0(r) = 0b + Q$ c - 3^ (0<r<r c ), (44) 

0(r) ^ 4> B --lQ~^-e- mB ^- r ^ (r>r c ). (45) 
r 

From Eq. the effective coupling outside of the body is of the order of Q, which means that local gravity constraints 
are not satisfied unless \Q\ is very much smaller than unity. 

B. Thin-shell solutions ((r c - ri) < r c ) 

From Eq. (jltl)) we find that it is possible to make the effective coupling small if n is close to r c . This corresponds 
to the thin-shell regime in which the field is stuck in the region < r < n with (r c — r\)jr c <C 1. Equation ((10)) 
shows that the mass term uia also affects the strength of the effective coupling. In the following we shall discuss the 
cases of large and small mass limits separately. 

1. The massive case (mAr c 3> 1) 

As we see below, thin-shell solutions originally derived in Refs. [l5|, [l6| can be recovered by taking the limit 

m A ri > 1 , (46) 

together with the thin-shell condition given by 

Ar c = r c - n < r c . (47) 
Expanding Eq. (|39[) in terms of small parameters Ar c /r c and l/m^fc, we obtain 

£th ^ + ^. (48) 
r c raATc 

When TfiA^c ^ (Ar c /r c ) _1 , this recovers the relation e t h — Ar c /r c derived in Refs. p^l. Il6t. 
From Eq. (|40|) the field profile outside of the body is given by 
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Neglecting second-order terms with respect to Ar c /r c and l/m^r,-, we get 

<j>(r) ~ 4> B 

where Q e ff is the effective coupling given by 



Qc« — 3Q 



Ar r 



m A r c 



3Qeth ; 



(50) 



(51) 



where we used the relation (|48[) . As long as both Ar c /r c and l/(mAf c ) a re much smaller than unity so that eth <C 1, 
it is possible to satisfy local gravity constraints. 

On using Eqs. (j3"5)) , (|36[) and (j4"i?|) together with Eq. we find that the field value at r = 0, r = r%, r = r c are 
given, respectively, by 



0(0) ~ 0A 

0(ri) ~ 4> A + 



12Q$ C 



6Q$ C 

{m A r c ) 2 ' 



(r c ) ~ A + 6Q$ C 



(m A r c 



I Ar c 

m A r c r c 



1 /At, 

2 I r r 



(52) 
(53) 

(54) 



Due to the presence of the e mAr<! term in Eq. (|5"2"]) . 0(0) needs to be very close to 0a- Under the condition m A r c 3> 1, 
the field can rapidly roll down the potential unless it is very close to the potential maximum at = <p A . Note also 
that when Q < (including f(R) gravity) we have 0(0) < 4> A (< 0) from ([52]) . Hence the boundary condition at 
r = that corresponds to thin-shell solutions with m A r c 3> 1 is chosen so that the field does not reach the curvature 
singularity at = [32l ]. 

From Eqs. (|52")) - (f5"4"|) it is clear that |0| gradually increases for larger r. The field stays around = 4> A within the 
body and it changes to the value 0s for r > r c . 



2. The light mass case (iriAr c <§C 1) 



Let us next consider the light mass case satisfying the condition 

m A ri < 1 , 



(55) 



together with the thin-shell condition Ar c <C r c . Note that this case is not equivalent to the thick-shell case in 
which the limit t\ — > is taken in Eqs. (j35| )- (|37jl . Since the radius T\ is close to r c , the ri-dependent term such as 
Qp A r\j (3r) in Eq. (|36j) is non- vanishing. 
From Eq. (|3l?|) we get the following relation 



Arv 



eth 



+ 



(m A r c ) 2 {m A r c ) 2 



(56) 



which gives the relation eth 3> 1- Even if the body has a thin-shell (Ar c /r c -c 1), the parameter eth is much larger 
than unity under the condition m A r c <C 1. The field profile (|40p in the region r > r c reduces to 



2Q 



1 - -^(m A r c ) 2 



GAL 



j(r- 



(57) 



This means that the coupling is of the order of Q as in the thick-shell case. Hence it is not possible to be compatible 
with local gravity constraints for \Q\ = 0(1). The effective coupling Q cS is of the order of Q even if r x exists in the 
intermediate regime (0 < r\ < r c ) away from r c . 
The field value at r = r r is estimated as 



3 ) ~ S - 2Q$ C ~ 4> B 



(58) 



where we used the relation 4>b — 6Q < i ) c /(m J 4r c ) 2 that comes from Eq. (|56[) under the condition \4>b\ 3> \<Pa\- Unlike 
thin-shell solutions with a massive limit (m^r c 3> 1), the field stays around ~ 0^ during the transition from the 
inside to the outside of the body. This property is different from the thin-shell case with the massive limit in which 
the field changes from <\>a (0 < r < r c ) to cf>B {r > r c ). 
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C. No-shell solutions (n = r c ) 



Let us next consider the case in which the solution is described by Eq. (|35[) in the whole internal region of the body. 
In this case Eq. (|3"5)) is joined to Eq. ([57)1 at the surface of the body (r = r c ). Setting r % = r c in Eqs. ([55]) and (|37[) . 
we get the following field profile: 

(0 < r < r c ) , (59) 



(r>r c ), (60) 



e -mB(r—r c ) 



' m A {e- mAr " + e mAr ") 
which agrees with the result in Ref. [3l|. We require that the following condition is satisfied: 

m A [<X r c) - 0a] < Qpa , (61) 

which is equivalent to 

eth < % — v ■ (62) 

Notice the essential difference from the "zero-shell" limit of the thin-shell case where the equality holds in Eq. (|G2l) . 

1. The massive case (mAr c S> 1) 

If the field <f> is massive such that ttiaTc ^> 1, the solution outside of the body is 

4>{r) ~4> B - 6Q— e th (l - e- m ^ r - r ^ . (63) 

r \ m A r c J 

This shows that the effective coupling is given by (Jeff — 3Qeth, which is the same as Eq. (fSTj) in the thin-shell case 
with the massive limit {ijiaTc S> 1). It is possible to satisfy local gravity constraints provided that eth *C 1. The field 
values at r — and r = r c are 

«r„) * A, + . (65) 

Hence 0(0) needs to be very close to <f>A- 
Equation ([62]) gives the following constraint 



eth < — ■ (66) 

TTlATc 

From Eq. (|48|) wc find that the opposite inequality, eth > 1/(™a?"c), holds for the thin-shell case in the massive limit 
(m A r c > 1). 

2. The light mass case (jtiaTc <C 1) 

When the field is almost massless such that m.AT c <§C 1, the solution in the region r > r c is 

CM 

<j>(r) ~ B - 2Q c -e th (m A r c ) 2 e- m ^ r - r ^ . (67) 

r 

The field stays around the value 4>b both inside and outside of the body. 
From Eq. ([62]) we get 

1 

{mATcf 

As long as e t h is much smaller than l/(m^r c ) 2 , it is possible to make the effective coupling Q c ff = Q eth(wA^c) 2 small. 
However, as we will see in the Appendix, the mass uia is generally heavy to satisfy the condition eth(^A^c) 2 3> 1 in 
concrete models that satisfy local gravity constraints. Hence this case is out of our interest. 



eth < , , 2 • (68) 
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V. CONCRETE MODELS 



In this section we consider experimental bounds on model parameters in concrete scalar-field potentials and inves- 
tigate how these models are consistent with local gravity constraints. In the previous section we have shown that the 
amplitude of the effective coupling Q ff can be made much smaller than \Q\ in two cases: (i) thin-shell solutions with 
iTiAfc ^> 1 and (ii) no-shell solutions with m,Ar c ^> 1. In both cases we have found 

QcS ^ 3Qe t h = ^ A ' ^ 69 -' 

which shows that |Q e ff| "C 1 under the condition \cf>B — 4>a\ <C $ c . 

The presence of the fifth-force interaction mediated by the field </> gives rise to a modification to the spherically 
symmetric metric. On using the thin-shell solution (|50p . the post Newtonian parameter 7 is given by 7 ~ (1 + 
"\/6Qcff/3)/(l — \/6<3cff/3) The present solar-system constraint on 7 is I7 — 1| < 2.3 x 10 -5 [29j . which comes 

from a time-delay effect of the Cassini tracking. This translates into the condition 

4.7 x 10~ 6 , x 
eth,® < — j^p- , (70) 

where eth,© is the thin-shell parameter for the Sun. Under the condition \4>b \ \4>a\, the bound (f7U|) corresponds to 

|0 S , | <5.9x 10- 11 , (71) 

where we used the value "I> c ,© ~ 2.1 x 10~ 6 . 

The fifth force induced by the field 4>(r) leads to the acceleration of a point particle given by = \Q c s4>( r )\ 0- 
This then gives rise to a difference for free-fall accelerations of the moon (aMoon) and the Earth (a©) toward the Sun. 
Using the present experimental bound, 2|aMoon — a©|/(aMoon + a©) < 10~ 13 [29|, under the conditions that the Sun, 
Earth, and Moon are subject to the effect of reducing the effective coupling discussed above, we obtain the following 
constraint [T(| 

8.8 x 10" 7 

eth,©< , (72) 

where eth,® is the thin-shell parameter for the Earth [lj|. The thin-shell condition for the atmosphere of the Earth 
provides the same order of the upper bound, provided that e t h — Ar c /r c [2Cl j . Using the value $ Ci © ~ 7.0 x 10~ 10 for 
the Earth, the bound ([72")) can be expressed as 

|0b,©| < 3.7 x 10- 15 . (73) 

We recall that 4>b,® (<^b,®) depends upon the density far outside of the Earth (the Sun). For both the Earth and 
the Sun we take the value pb — 10 -24 g/cm 3 that corresponds to the dark matter/baryon density in our galaxy. Then 
|<Ab,©| is equivalent to |</>b,©| provided that the scalar-field potential is specified. In the following we adopt the severe 
bound ([75)1 for two potentials given in Eqs. ([5]) and (Q2 



A. Inverse-power law potential 

Let us consider the inverse power- law potential ©. In this case we have 



n M p 4 / M x " +4 



1+1 



M pl . (74) 

Q p B \M p i 1 

On using the bound (|73p with n and Q of the order of unity, we get the following constraint 

M< lQ-^^-Mpi. (75) 

This shows that M < 10" 2 eV for n = 1 and M < 10" 4 eV for n = 2, which are consistent with the bound derived in 
Ref. 0. 
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The mass squared m A about the potential minimum at <f> = 4>a is given by 

m A = n(n + l)\ 




Multiplying the squared of the radius r c of the Earth and eliminating M with the use of Eq. (|74|) . we obtain 

i 

(m A r c ) 2 =6(n + l)Q$ c , e [ ^ ) + ' — !— . (77) 

\PB J <PB,® 

Since the mean density of the Earth is p A ~ 10g/cm 3 , the experimental bound ([75)) leads to 

m A r c > 3yJ(n + l)Q ■ . (78) 

When Q is of the order of unity one has m A r c ^ 10 9 for n = 1 and 771^4?^ > 10 7 for 77 = 2. Hence in these cases the 
mass niA in fact satisfies the condition m A r c ^ 1 so that the effective coupling is given by Eq. (|69|) . 

If the body has a thin-shell then it is required that the condition. e t h > 1/(777,47^), is satisfied. Employing the 
relation ([77|) . this condition corresponds to 

m A r c < (n + 1) ( — ) " + " ~ (n + 1) • 10^ . (79) 



For the compatibility of this inequality with Eq. (|78p . we require that 

n<4, (80) 

in which case the body has a thin-shell. Meanwhile, when n > 4, the solution corresponds to the no-shell case with 
mAr c > 1. 

We note that the condition ([79")) for the existence of thin-shells can be also expressed as 

6Q4> e , ffl / / Pb\^ t rjm 

0s,© > — — ■ (81) 

n + 1 V / 

When Q = 0(1) we have 4>b,® ^ 10~ 21 for n = 1 and </>b,© > 10 -17 for n = 2. If future experiments provide tight 
bounds on </>b,© such as <^b,© ^ 10~ 21 , the field profile corresponds to the no-shell case for n > 1. Hence the body 
does not have a thin-shell in such a situation. We stress again that no-shell solutions are also consistent with local 
gravity constraints, since the effective coupling Q c s takes the same asymptotic form as that in the thin-shell case with 
the massive limit. 



B. The potential motivated by f(R) gravity 

The next example is the potential (Tl2|) . which covers viable f(R) models that satisfy local gravity constraints 
[1, @, The field value 4>b,® is given by 

Note that Vq is of the order of the present cosmological energy density po ~ 10 _29 g/cm 3 , if the potential of the 
field <fr is responsible for a late- time acceleration [l(|. Employing the experimental constraint (|75|) with C and |Q| of 
the order of unity, one can derive the bound p > 14/9. More precise analysis using the information of the late-time 
acceleration provides the constraint [l(| : 

P>1 ~13.8-l 5 og 10 |Qr ^ 



When \Q\ = 0.1 and \Q\ = 1 this gives the bounds p > 0.66 and p > 0.64, respectively. 
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The quantity (mAT c ) 



2 for the Earth is 



(m A r c ) 2 =6|Q|(1 



P) 



|0B,©| \PB J 



(84) 



Under the constraint |73|) we get 



nw c >[|Q|(l 



-I / 31— 6p 



(85) 



When p > 0.65 this condition correspond to ttlaTc ^ 10 39 , which means that the field is extremely massive inside of 
the body. 

The necessary condition under which the body has a thin-shell, eth > 1/(™a^c), is given by 



It is clear that the values of tuaTc that satisfy both (|85|) and (|86|) exist for p ranging in the region < p < 1. Hence 
the body has a thin-shell under the current experimental bound (|73p without reaching the no-shell regime. 
The condition (f8"6")) can be also expressed as 



When p = 0.65 and \Q\ — 1 one has |</>,b,©| > 10 -79 . Note that this field value gets even smaller for p larger than 0.65. 
Since it is unlikely that future experiments reach the level of the constraint < 10 -79 , the field profile is well 

described by thin-shell solutions for realistic bounds on |^>_b,©|- This situation is different from that for the inverse 
power-law potential. 



In this paper we have derived analytic solutions of a chameleon scalar field </> in the background of a spherically 
symmetric body by taking into account a field mass about the potential minimum at <fi = cf>A inside of the body. 
The mass ttia is important to determine the field value <fi(ri) at which <p begins to evolve along the potential. This 
is characterized by the condition, m 2 A [4>{ri) — 4>a] — Qpa, where Q is a strength of the matter coupling and pa is a 
mean density of the body. In the region r\ < r < r c (r c is a radius of the body) , the presence of the matter coupling 
leads to a rapid evolution of the field. It is known that the effective matter coupling outside of the body gets much 
smaller than \Q\ if the body has a thin-shell (l5l. [l6|. 

By considering the solution in the region < r < r\ with the mass tua taken into account, we have derived Eqs. (135D - 
(|37[) as analytic expressions of the field profile. There exist three different regimes depending on the position of ri, 
i.e., (i) thick-shells (r\ = 0), (ii) thin-shells ((r c — r\)jr c <C 1), and (iii) no-shells (r\ = r c ). In the thick-shell case we 
have recovered the usual solution for <f> in the exterior of the body, where the matter coupling is of the order of Q. 

The standard thin-shell solution derived in Ref. [TBI fl6| is recovered in the limit thaTc — * oo. When toa^c S> 1 
the effective coupling outside of the body is given by Q c g ~ 3Q(Ar c /r c + l/mAr c ) and is thus subject to the 
correction coming from l/mAr c . However, if one uses the thin-shell parameter e t h = {4>b — 4>a)I{QQ^c) introduced 
in Refs. 0,0], this effective coupling can be expressed by Q e ff — 3Qe t h- Since the condition e t h -C 1 is used when 
we place experimental bounds on model parameters, the original works [l5l.[l6| based on this criterion are not subject 
to change by the inverse mass correction. When m,Ar c <C 1 the effective coupling in the exterior of the body is of 
the order of Q, which means that models with large couplings (|Q| = 0(1)) are not consistent with local gravity 
constraints. The requirement of the large mass (mAT c 3> 1) for the consistency with local gravity experiments implies 
that the field needs to change rapidly from 4>a to 4>b during the transition from the thin-shell regime to the exterior 
region of the body. 

In the no-shell case, the mass term m^(<^(r c ) — <Pa) is smaller than the matter coupling term Qpa so that the field 
always stays around cf> = 4>a inside of the body. The interior and exterior solutions of the field arc given in Eqs. (|59p 
and (|6TJ|) . We have shown that, in the limit toa^c 3> 1, the field profile outside of the body and the expression of the 
effective coupling Q c g — 3Qeth are the same as those in the thin-shell case with the same massive limit. Since there is 
no lower bound for e t h contrary to the thin-shell case, there is a possibility that the no-shell case is compatible with 
local gravity constraints even in the situation that the thin-shell case is not. 

We used experimental bounds coming from the violation of equivalence principle as well as solar-system tests to 
constrain concrete models under the chameleon mechanism. For the potential V(4>) = M 4+n (f>~ n the experimental 



m-ATc < (1 — p) ■ 10 1 ~ p ■ 



(86) 




(87) 



VI. CONCLUSIONS 
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bound (|72|) leads to the constraint m,AT c ^> 1 for the Earth (e.g., vtiaTc ^ 10 9 for n = 1 and Q = 1). Under the current 
experimental bound we have found that the body has a thin-shell for n < 4, whereas the field profile corresponds to 

_ 15re+130 u 

no-shell solutions for n > 4. As long as the condition, M < 10 ™+ 4 pl , is satisfied, the inverse power-law potential 
is consistent with local gravity tests both for thin-shell and no-shell solutions. 

We also considered another potential given in (fT2"|) that includes viable f{R) models as specific cases. In this case 
the field is very massive inside of the body: mAr c > 10 39 for the values of p consistent with local gravity constraints 
(P ^ 0.65). We have found that under the current experimental bound the field profiles always correspond to thin-shell 
solutions instead of no-shell solutions. This different property compared to the inverse-power law potential comes 
from the fact that the field inside of the body is much heavier. 

We have thus shown that the chameleon mechanism works in a robust way provided that the field mass inside of the 
body satisfies the condition ttlaTc S> 1 with appropriate boundary conditions. In both thin-shell and no-shell cases 
we just need to use the effective coupling Q c r = 3Qeth = {4>b — 4>a)/^c to place experimental bounds on model 
parameters. We note that under the condition ttiaTc 3> 1 the field value at r = is required to be very close to 4>a, 
see Eqs. (|5"2l and ([M]) . This property is especially severe for the potential fT2]) due to the condition niAr c > 10 39 . 

It will be interesting to extend our analysis to the case in the strong gravitational background as in the recent work 
[HI of f(R) gravity (motivated by Ref. [32} )■ Under a strong gravity we require a careful analysis by taking into 
account the effect of gravitational potentials to the field equation in addition to the issue of the boundary condition 
at r = 0, which we leave for future work. 
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Appendix 

In this Appendix we show that no-shell solutions with a nearly massless field does not satisfy experimental bounds 
for the models discussed in Sec. El In this case the effective coupling is given by Q e g = <5eth(m- J 4' , c) 2 from Eq. (|oT)) . 
The experimental test for the violation of equivalence principle using the free-fall acceleration of the Moon and the 
Earth toward the Sun provides the bound |Q c ff | < 2.6 x 10~ 6 . This translates into the condition 

\^ B ,®\{m A r c ) 2 < 1.6 x 10- 5 $ c ,© . (88) 

For the inverse power-law potential ([5J, the constraint ((55)) leads to 

Q{n+ 1) (jf^ n+1 < 2 - 7 x 10 ~ 6 ' ( 89 ) 

where we used Eq. (|77j) . Since pa/ Pb — 10 25 for the Earth, it is not possible to satisfy the inequality (|89)) for n > 
and Q = 0(1). 

For the potential (fT2|) . the constraint ([88]) leads to 

IQI(I-P) (^)" P <2.7x 10- 6 , (90) 
where we used Eq. ((54"|) . Again we can not satisfy this condition for < p < 1 and \Q\ = 0(1). 
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